Excitons in Molecular Aggregates with Levy Disorder: 
Anomalous Localization and Exchange Broadening of Optical Spectra 
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We predict the existence of exchange broadening of optical lineshapes in disordered molecular 
aggregates and a nonuniversal disorder scaling of the localization characteristics of the collective 
electronic excitations (excitons). These phenomena occur for heavy-tailed Levy disorder distribu- 
tions with divergent second moments - distributions that play a role in many branches of physics. 
Our results sharply contrast with aggregate models commonly analyzed, where the second moment 
is finite. They bear a relevance for other types of collective excitations as well. 
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During the past decade much work has been devoted 
to low-dimensional optical materials, such as molecular 
aggregates, photosynthetic antenna systems, conjugated 
polymers, and quantum wells and wires. They derive 
their interest from strong and color-tunable absorption 
and luminescence properties, fast electromagnetic energy 
transport, or optical switching Common to this wide 
variety of materials is that their properties are rooted in 
excitonic eigenstates, collective excitations that consist 
of a linear combination of local excited states. The low- 
dimensional structure, however, makes these excitations 
strongly dependent on the presence of disorder in the 
site energies and (or) the intersite interactions. Disorder 
localizes the excitons and thereby affects basic proper- 
ties, such as the optical oscillator strength per state, the 
absorption linewidth, and the excitation diffusion con- 
stant. The standard disorder model considers Gaussian 
or box-like distributions for the energies or the interac- 
tions, which both have a finite second moment. 

In this Letter, we show that disorder models with di- 
vergent second moments give rise to drastically different 
results for the exciton states and the collective optical 
response than the standard ones. To this end, we study 
molecular aggregates with site disorder taken from distri- 
butions that have caught interest in a wide range of phys- 
ical problems, namely heavy-tailed Levy distributions |2j. 
These were introduced originally as "exceptional cases" 
that do not obey the central limit theorem. However, 
during the past 20 years it has been recognized that they 
frequently occur in physics and their divergent second 
moments give rise to a variety of surprising effects in sub- 
fields ranging from statistical physics to optics, plasma 
physics, and condensed matter physics |3l4l3|. In the con- 
text of the model studied in this Letter, it is of particular 
interest that molecules embedded in a structurally disor- 
dered (glassy) host provide an interesting realization of 
systems with heavy-tailed disorder. The interactions be- 
tween the molecule and the multipoles representing the 
structural disorder of the host give rise to random shifts 



of the energy levels of the molecule, which is reflected 
in the broadening of the absorption line. Indeed, it has 
been shown experimentally and theoretically that these 
absorption lines have heavy tails (decaying slower than 
Lorentzian) (Tol EL with mean and variances obeying 
Levy distributions |l2l. [l3| . 

In previous studies, Levy statistics applied to a single 
random variable, such as the displacement of a diffusing 
particle or excitation energy of a single molecule. Here we 
present the first study of interacting degrees of freedom 
that individually obey heavy-tailed Levy statistics and 
predict totally new phenomena. More specifically, we 
show that the model considered has remarkable collective 
optical properties that differ even qualitatively from the 
well-studied case of Gaussian disorder 14]. In particular, 
we find the phenomenon of exchange broadening of the 
absorption line shape, the counterintuitive appearance 
of fine structure in the density of states (DOS) and the 
absorption spectrum with increasing disorder strength, 
and a nonuniversal disorder scaling of the distribution of 
exciton localization lengths. 

We consider a one-dimensional Frenkel exciton model 
to describe the optical properties of linear aggregates of 
dye molecules. In fact, this model is very generic and also 
is used often to describe the properties of the other opti- 
cal materials mentioned in the introduction. The model 
consists of a linear array of two-level molecules with par- 
allel transition dipoles, which interact through dipole- 
dipole transfer interactions. The optical excitations of 
this system are described by the eigenstates of the exciton 
Hamiltonian matrix H nm = E n 5 nm - J (5 m ,n+i + S m , n -i) 
(l4| , where for simplicity only nearest- neighbor couplings 
are considered [15|- The transfer integral J > is as- 
sumed to be constant. Disorder is included by taking 
the site energies E n as uncorrelated stochastic variables, 
drawn from a distribution p(E). 

Special for our model, as compared to previous studies, 
is that we consider heavy-tailed Levy distributions for 
p(E). In particular, we consider the class of symmetric 
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Levy distributions with mean zero, given by 16 1 

1 r°° 

P( E )=2^J dte lEt exp(~\atn . (1) 

Here, a > and < a < 2 are called the scale pa- 
rameter and index of stability, respectively. The former 
determines the half-width at half maximum (HWHM) 
of p(E) and we will hereafter refer to it as the disorder 
strength, while the latter fixes the asymptotic behavior 
p(E) ~ l/E 1+a for E > a and a < 2. The power-law 
behavior for large E gives rise to heavy tails in p(E), 
which in turn lead to a divergent second moment. 

It should be noted that the average site energy for a 
disorder realization on the chain, E = N^ 1 Yl n =i 
also obeys a Levy distribution, with the same index of 
stability, but a renormalized disorder strength [ltjj ]: 



a* = aN- 



(2) 



This property is referred to as stability of the distribu- 
tion. For a Gaussian distribution (a — 2), Eq. ((2]) gives 
the well-known result a* = a/VN, which reflects the 
exchange narrowing effect of the optical line shape of 
molecular J aggregates [17] : the energy distribution of 
delocalized excitons is narrower than a, because they av- 
erage over the energy fluctuations of the individual sites. 
By analogy one expects from Eq. ([2]) that for a < 1, this 
effect is absent. The Lorentzian case (a = 1) is special, 
as it yields a* — a. This case was studied in Refs. [lij : 
while it does not exhibit exchange narrowing, its main 
properties still are similar to those for Gaussian disor- 
der. Here, we are particularly interested in distributions 
with a < 1, for which a* > a (exchange broadening). 

We analyzed our model by numerical simulations. We 
used the algorithm described in Refs. to generate the 
random energies E n , after which the Hamiltonian was 
diagonalized numerically and the exciton DOS, the ab- 
sorption spectrum, and the exciton localization lengths 
N\ oc (inverse participation ratio) were calculated using 
standard methods [14j. We used chains of 200 sites and 
averaged over tens of thousands of disorder realizations. 
As typical examples, the results obtained for a = 1/2 are 
presented in Figs. Q]- |U 

Figure [1] shows the DOS and the absorption spectrum 
for various disorder strengths a. The DOS exhibits a typ- 
ical one-dimensional (ID) shape, with peaks at the band 
edges E — ±2 J that are smeared by disorder. However, 
three additional features are visible: one at the band 
center, E = 0, and two more at E = ±J. The relative 
importance of these three features strongly depends on 
a and a, revealing a transition from a ID excitonic DOS 
to a mostly monomeric DOS with increasing a. The ori- 
gin of the extra features, which we found to get more 
pronounced upon decreasing a, will be explained below. 

The structure in the DOS is also reflected in the ab- 
sorption spectrum [Fig. [lib)]. While this spectrum is 
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FIG. 1: DOS (a) and absorption spectra (b) for Levy disor- 
der with a = 1/2 and disorder strengths a = 0.01J, 0.1J, 
0.2J, 0.5J, and J. The arrows mark the curves in order of 
increasing a. 



dominated by the intense band edge peak (the J band), 
characteristic for J aggregates, two much less intense fea- 
tures occur at E = — J and E = 0. For Gaussian or 
Lorentzian disorder these features cannot be discerned. 
Figure[ljb) also reveals another new effect: with increas- 
ing value of a, the J band shifts to the blue, while for 
Gaussian disorder a red-shift occurs 
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FIG. 2: The HWHM of the absorption spectrum as a function 
of a for Levy disorder with a — 1/2 (crosses) and for Gaus- 
sian disorder (circles). The corresponding power law fits are 
plotted as a solid and a dashed line, respectively. The thin 
solid line HWHM = a is plotted for reference. 

Figure [5] gives the HWHM of the absorption band as a 
function of a, for a = 1/2 (crosses) and for comparison 
also for a — 2 (Gaussian; circles). This plot confirms the 
most fascinating effect of heavy-tailed Levy distributions, 
anticipated above, namely the occurrence of exchange 
broadening: While for Gaussian disorder the HWHM is 
smaller than the bare disorder er, for a = 1/2 it is larger 
(for < a < 0.6). The best power-law fit to the data 
reads HWHM = 0.85 J (cr/J)°' 60±0 '° 3 (solid line), which 
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differs strongly from the scaling HWHM oc J (a j J) 4 ^ 3 for 
the Gaussian case (dashed line) [14j - While the broad- 
ening effect qualitatively agrees with Eq. ([2]) , we will see 
below that the quantitative explanation of the exponent 
0.60 ± 0.03 is more subtle. 



to a* given in Eq. © with N replaced by N* [HJ]. Gen- 
eralizing to the case of Levy distributions yields 




FIG. 3: The average, iVioc, (crosses) and the standard de- 
viation, 5N\ oc , (dots) of the localization length distribution 
P(JVi oc ) as functions of a for Levy disorder with a = 1/2. 
The thick solid and dashed lines are power-law fits, respec- 
tively (see text). The thin solid line gives 10(5N\ OC /Ni oc ), 
while the inset shows the localization length distributions for 
a = 0.001J (dashed) and for a = 0.1J (solid). 

Finally, Fig. [3] characterizes the distribution P(N\ oc ) 
of the localization length of the exciton states that oc- 
cur in the energy interval [— 2.1J, — 1.9 J], i.e., around 
the lower band edge Et, = — 2J, the region that domi- 
nates the optical response. Shown are the numerical data 
(symbols) and their power-law fits (lines) for the average, 
N\ oc , (crosses and solid line) and the standard deviation, 
SNioc, (dots and dashed line) of P(iVi oc ) as a function 
of a. The fits represent the data very well, and read 
Moc = 2.56(J/ct)_ - 36 and 8N ioc = 0.91 (J/ct) 43 . Thus, 
the ratio 5N\ ac / N\ oc is not constant, implying that, in 
contrast to Gaussian and Lorentzian disorder [18| , Levy 
disorder with a < 1 does not lead to a universal function 
for P(N\ oc ) in the cr-range considered. This is clearly ob- 
served from the inset of Fig. [3l where P(Ni oc ) is plotted 
for a = 0.001J and a = 0.1J, after scaling the horizon- 
tal axis to the average localization length appropriate for 
the a value considered. We note that the scaling rela- 
tions described here, only weakly depend on the choice 
of the energy interval around E^. 

The explanation of all phenomena found above, lies 
in the interplay between two mechanisms for localizing 
the exciton states around Eh by heavy-tailed Levy dis- 
order. For Gaussian disorder only one such mechanism 
exists, namely localization of states in effective potential 
wells created by random site energies [2(|. The typical 
localization length TV* of these states can be obtained 
from a balance of the energy spacing ~ J/N* 2 of two 
adjacent states within a localization segment and their 
disorder-induced scattering rate, which scales according 



N* = 



( Sir 2 J 



(3) 



where £ Q is a numerical factor of order unity (£1/2 ~ 0.7). 
For a = 1/2, this gives N* cx {J/a) 1/3 . Generally, N* 
is not identical to 7Vi oc because the latter also includes 
effects of segmentation (see below). 

For Gaussian disorder, N* is the only length scale rel- 
evant to the band edge states and the optical response. 
In the case of the heavy-tailed Levy distributions, how- 
ever, a second localization mechanism - and correspond- 
ing length scale - exists. The long tails lead to a high con- 
centration of outliers, i.e., sites with energy \E n \ > 2 J. 
These fluctuations are so large that the interaction J can- 
not overcome them; they therefore break up the chain in 
segments of length N scgl capped by two outliers, which 
form the maximum intervals over which excitons may 
delocalize. For decreasing values of a this effect gets 
stronger. It is straightforward to show that for a < 1 the 
segment length distribution is exponential, with mean 



2T(a) sin(7ra/2) 
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(4) 



iV S eg is the second length scale in the problem. 

The existence of two localization mechanisms is con- 
firmed by Fig.|H which for a — 1/2 and a = 0.1J shows a 
typical realization of the exciton wave functions and ener- 
gies in the neighborhood of the lower band edge and well 
below it. Indeed, we observe quite a few outliers, with 
positions indicated by the vertical dashed lines. Each 
outlier has a strongly localized s-like, i.e., without nodes, 
exciton state (lower panel); not all these states are seen, 
as several of them have energies outside the plot's range. 
With regards to the two localization mechanisms, three 
situations may be distinguished, each of which occurs in 

Fig.m 

(i) TV* < iV se g. This is the situation common for Gaus- 
sian disorder. Near the lower band edge one then typi- 
cally finds doublets of s- and p-type states (with no nodes 
and one node, respectively) roughly localized on the same 
interval [2jJ. In Fig. H this occurs between n — 52 and 
n = 96. 

(ii) N* pa A^cg. Here one finds multiplets of three 
or more states on a single segment, which resemble the 
states of a disorder- free chain of size iV SG g. In Fig. 0] 
examples occur between n — 1 and n — 31, between 
n = 92 and n = 112, and between n = 164 and n = 200. 

(iii) N* > N scg . In this case the segmentation strongly 
confines the excitons and the states within a segment 
typically get further separated than in the absence of 
segmentation. If the segments are still relatively large, 
their lowest s-like eigenstates contain considerable oscil- 
lator strength and occur just above the band edge (e.g., 
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FIG. 4: Exciton wave functions and energies for a typical 
realization of Levy disorder with a — 0.1J and a = 1/2. 
The upper panel focuses on energies around the lower band 
edge; the lower panel displays an energy interval deep in the 
DOS tail. The dashed vertical lines indicate the positions of 
outliers in the site energy. 

the state between n = 154 and n — 164 in Fig. 2]). With 
increasing value of a, the segments get shorter and the en- 
ergy of the s-like state grows. This explains the blue-shift 
of the J band found in Fig.QJb). Also, for growing a, seg- 
ments of length 1 and 2 become more likely. These give 
rise to states with energies distributed around the average 
monomer and dimer energies, E = and E = ±J, re- 
spectively, thus explaining the extra features in the DOS 
and the absorption spectra observed in Fig. [1] Closer 
scrutiny even reveals features for segments of length 3, 4, 
etc, but these are weak and overshadowed by the band 
edge peaks. As is seen in Fig. [TJ all these extra fea- 
tures are washed out when a approaches J, because they 
broaden and cannot be distinguished anymore. 

Of course, the existence of two localization mecha- 
nisms, each with its own length scale, is directly re- 
sponsible for the observed nonuniversality of the local- 
ization length distribution. Moreover, the interplay be- 
tween both mechanisms also explains the disorder scal- 
ing of the exchange broadening of the absorption band. 
Using Eq. ^ and replacing N by either N* or iV SC g, 
we obtain the contributions to the HWHM from the 
states localized by the two different mechanisms, respec- 
tively. This leads to HWHM* cx J(a/ J) 2a ^ 1+a ^ and 
HWHM SC g cx J{a/J) a . For a = 1/2, the scaling expo- 
nent equals 2/3 and 1/2, respectively. Since the numeri- 
cally obtained exponent equals 0.60 ± 0.03, we conclude 
that both scales, N* and iVseg, almost equally contribute 



to the disorder scaling of the J band width. 

In summary, we have shown that heavy-tailed Levy 
disorder has dramatic consequences for the optical prop- 
erties of linear molecular aggregates. Novel effects occur, 
such as exchange broadening and a nonuniversal scaling 
of the distribution of exciton localization lengths. We 
have shown that these effects may all be traced back to 
the simultaneous occurrence of two different localization 
mechanisms, each with its own length scale. We expect 
that the occurrence of two localization length scales has 
equally dramatic effects on other collective excitations 
subjected to heavy-tailed disorder. 
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